Explanation of Di-jet asymmetry in Pb+Pb collisions at the Large Hadron Collider 
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We investigate the medium modification of a partonic jet shower traversing in a hot quark-gluon 
plasma. We derive and solve a differential equation that governs the evolution of the radiated gluon 
distribution as the jet propagates through the medium. Energy contained inside the jet cone is 
lost by dissipation through elastic collisions with the medium and by scattering of shower partons 
to larger angles. We find that the jet energy loss at early times is dominated by medium effects 
on the vacuum radiation, and by medium-induced radiation effects at late times. We compare our 
numerical results for the nuclear modification of the di-jet asymmetry with that recently reported 
by the ATLAS Collaboration. 



Quarks and gluons scattered to large transverse mo- 
mentum pt during the early stage of a relativistic heavy- 
ion collision are regarded as powerful probes of the tran- 
sient quark-gluon plasma produced in these reactions [J . 
Such hard partons interact with the traversed matter and 
lose some of the initial energy by elastic and inelastic col- 
lisions with constituents of the medium 0-01 • The essen- 
tial aspects of parton energy loss have been confirmed by 
observations of strong high-p^ single hadron suppression 
in nucleus-nucleus collisions at the Relativistic Heavy- 
Ion Collider (RHIC) [E, 0- Other observations support- 
ing this interpretation include the suppression of back- 
to-back di-hadron and correlated back-to-back photon- 
hadron emission 0-0 ■ These phenomena are commonly 
referred to as "jet quenching". Although direct calori- 
metric measurement of jets has been performed at RHIC, 
detailed studies of the medium modification of jet show- 
ers are constrained by limited center-of-mass (cm.) en- 
ergy [iMl. 

The launch of the heavy-ion program at the CERN 
Large Hadron Collider (LHC) has extended the available 
cm. energy range by more than a factor of ten, enabling 
the investigation of medium effects on jets with trans- 
verse energies Et in excess of 100 GeV. Very recently, 
the ATLAS Collaboration reported first results on the 
modification of di-jet properties in Pb+Pb collisions at 
2.76 TeV 



^/snn ~ 2.76 TeV [13|. The most spectacular aspect 
of these results is the occurrence of extreme energy im- 
balances between two correlated jets emitted in opposite 
azimuthal directions around the beam axis. Similar pre- 
liminary results were reported by the CMS Collaboration 
[3] ■ These results provide a challenge for theoretical de- 
scriptions of the parton shower evolution in hot QCD 
matter. While a qualitative analysis of the physics pro- 
cesses responsible for the observed large di-jet asymmetry 
was presented in Ref . [T5| , no quantitative description of 
the experimental results has been reported to date. 

In this work, we study the medium modification of 
a parton shower in an expanding quark-gluon plasma 
and calculate the energy loss experienced by the shower 
within a cone ang le defined as i? = y/{A<j)y + {Arj)^. 
The total loss of energy from the jet cone is the sum 
of energy loss incurred by the leading parton and its 
radiated gluons, as well as the energy carried by glu- 



ons that are scattered out of the jet cone. The leading 
parton may transfer energy to the medium by collisions 
with matter constituents; this energy loss is quantified 
by the transport coefficient e — dE/dt, denoting the col- 
lisional energy loss per unit time. The leading parton 
also radiates gluons due to scattering in the medium; 
this effect contributes to the depletion of energy con- 
tained within the jet cone if the radiated gluons are emit- 
ted outside the cone. In addition, the radiated gluons 
(from both vacuum and medium-induced radiation) orig- 
inally contained inside the jet cone not only suffer colli- 
sional energy loss, but they may be kicked out of the jet 
cone by multiple scatterings in the medium after their 
emission. The transverse momentum broadening of the 
jet shower is characterized by the transport coefficient 
q ~ d{Ap±)'^ /dt [1, [l^. Here we relate these two trans- 
port coefficients hy q = 4Te, assuming a medium close 
to local thermal equilibrium and the applicability of the 
fluctuation-dissipation theorem [l3| • 

We now derive a differential equation that describes 
how the radiated gluon distribution /g(a;,fc^,t) = 
dNg{u},kj_,t)/dujdk'^ changes when the jet shower prop- 
agates through the medium. Here w denotes the gluon 
energy and k± the transverse momentum with respect to 
the jet axis. The resulting modified gluon distribution is 
used to calculate the energy loss from radiated gluons in- 
side the jet cone. The total energy loss from the jet cone 
is the sum of the energy loss from the leading parton and 
the radiated gluons, AE = AE'l + AEg. 

Consider a partonic shower that is propagating 
through a hot quark-gluon plasma. During their prop- 
agation, the radiated gluons will transfer energy into the 
medium by elastic collisions and accumulate transverse 
momentum in the process. The interaction between the 
leading parton and the medium may induce additional 
gluon radiation given by the rate dN™'^'^ / dujdk^dt. After 
their emission, these gluons may also lose energy and ac- 
cumulate transverse momentum. Combining these three 
contributions, we obtain the following differential equa- 
tion for the radiated double differential gluon distribu- 
tion: 



dt 



fg{uj,k\,t) = e 
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dN^""" 
dojdk^dt 
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Note that the uj dependence of the jet shower evolution 
has been discussed in Ref. [ll], where the medium re- 
sponse to the energy deposition of a jet shower was stud- 
ied (see also Ref. p^). 

The above equation must be supplemented with an 
equation for the initial gluon distribution, /g(w,fc^,ti) 
generated by vacuum showering of the primary parton 
before interacting with the thermal medium at time 
ti. The vacuum radiation spectrum is generated from 
PYTHIA H, and the gluon distribution /™'=(a;, fc^ , t,) 
at the initial time is obtained by requiring that only glu- 
ons with a formation time Tf = 2Ex{l — x)/k\ smaller 
than ti are radiated, where x = uj/ E denotes the energy 
fraction of the radiated gluons. 

In the higher- twist formalism of jet quenching, the rate 
for medium-induced gluon radiation is obtained as [2ll . 



dN^^^ 2as xP{x)q{t) . ^t-U 



diodkj_dt TT 



(2) 



where P(x) is the vacuum splitting function. The trans- 
port coefficient q depends on time because of the ex- 
pansion of the medium as well as the changing position 
of the primary parton. When solving for the radiated 
gluon distribution fgluj, k\,t), we impose a lower cut-off 
Wmin = 2 GeV on the radiated gluon energy, and any ra- 
diated gluon whose energy drops below Wmin is considered 
as part of the medium. This energy cut-off is about four 
times the highest temperature of the medium created in 
Pb-I-Pb colhsions at the LHC. 

As we already discussed, our initial condition for the 
radiated gluon distribution is taken as the vacuum radia- 
tion shower obtained at time U , after which we start the 
medium evolution of the gluon distribution. This is moti- 
vated by the observation that the vacuum contribution to 
the primary parton virtuality dominates at earlier time 
while the virtuality contribution from the medium takes 
over at later times |23| . Furthermore, interference effects 
among subsequent gluon emissions ensure that the radi- 
ated gluons in the vacuum are angular ordered, i.e. that 
gluons emitted at late times are radiated into small an- 
gles within the jet cone. Vacuum radiation at late times 
will thus not be sufficiently modified by the medium to 
contribute to the energy loss of the shower outside the 
jet cone. We note, however, that the vacuum radiation 
would need to be included as a source term in Eq. ([1]) 
in order to account for the interference between vacuum 
and medium-induced radiation. 

After solving Eq. ([T]) and obtaining the time evolution 
of radiated gluon distribution fg{to, k\,t), we may calcu- 
late the final energy of the gluons contained inside the 
jet cone: 



Eg{tf,R)= / LjdLjdklfg{uj,kl,tf), 



(3) 



where the subscript R indicates that the integration is 
performed over the interior of the jet cone with the con- 
dition k±/uj < R. The total energy inside the jet cone 



is the sum of the energies of the leading parton and the 
gluons inside the cone: 

E;,,{tf,R)^E^{tf)+Eg{tf,R), (4) 

where the final energy for the leading parton is given by 



Ei^itf) = Ei^iU)- J e{t)dt 
Lodujdk I dt 



diddk^dt 



(5) 



The above equation describes the time evolution of lead- 
ing parton energy in the medium. The momentum broad- 
ening of the leading parton in the medium may be in- 
cluded in a similar way. The leading parton energy at 
the starting time U is obtained as 



Ei^iU) = Ej,t{0) - / ujdujdkifg{uj,ki,U). 



(6) 



The complete reduction of energy contained within the 
jet cone is 



IS.E — Ejct {ti , R) — Ejct (tf , R) . 



(7) 



The initial energy within the jet cone is the sum of leading 
parton energy and the energy of the gluons inside the 
cone: 

E^^t{U,R)^EL{U)+ [ ujdujdklfg{uj,kl,U). (8) 

JR 

The final expression for the energy lost from the jet cone 
is: 

AE ^ [ LjdLjdkl[fgiu;,kl,t,)~ fg{uj,kl,tf)] 

JR 



ujdujdk , dt — 



dcudk^dt 



■f 



e{t)dt. (9) 



The first three terms describe the energy loss due to radi- 
ated gluons, while the last term is the coUisional energy 
loss by the leading parton. Equation (jH]) is quite general 
and can be understood if one splits the medium-induced 
radiation term into two parts: the radiated energy inside 
and outside the jet cone. Then the first and third term 
represent the gluon energy gains inside the cone, the sec- 
ond term is the remaining gluon energy inside the cone, 
and the last two terms represent the energy loss from the 
leading parton. 

We first apply the above formalism to the energy loss 
of a jet in a static medium. The result is shown in 
Fig. [1] where we show the contributions from vacuum 
radiated gluons (left) and medium-induced radiated glu- 
ons (right) to the total energy loss of the jet defined 
by a cone angle R = 0.4. We begin with a quark jet 
with initial energy E = 100 GeV which radiates glu- 
ons in vacuum before ti = \ fm/c, when we turn on the 
medium modification for the vacuum shower. The figure 
shows the additional energy loss experienced by the jet 
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FIG. 1: (Color online) The energy loss experienced by a 
100 GeV quark jet traversing a brick of quark- gluon plasma 
with constant temperature T = 250 MeV. Contributions due 
to different mechanisms of energy loss from vacuum radiation 
(left) and medium-induced radiation (right) are compared. 



shower inside the medium. Here we set the temperature 
of the medium to be T = 250 MeV and the transport 
coefBcient to be q = 0.7 GeV^/fm. The figure shows 
that vacuum radiation dominates the energy loss at early 
times; medium-induced radiation later takes over due 
to its much stronger length dependence. For both vac- 
uum and medium-induced radiation, we observe that the 
most significant contribution originates from collisional 
energy loss experienced by radiated gluons. For medium- 
induced radiation, the transverse momentum broadening 
of shower gluons presents another sizable contribution, 
while the contribution from radiation outside the jet cone 
is small. 

We now apply our formalism to calculate the jet cone 
energy loss and the medium modification of the di-jet 
asymmetry in Pb-|-Pb collisions at the LHC reported by 
the ATLAS Collaboration. To proceed, we scale the 
transport coefficients according to the temperature of 
the medium, q (x T^, with a constant factor adjusted 
to fit the experimental data. We model the space-time 
profile of the temperature T{f±,t) at midrapidity as fol- 
lows. The initial spatial profile for the entropy density 
s ^ T"^ of the medium is set to be proportional to the 
density of participating nucleons in the colliding nuclei. 
The nuclear density distributions are taken as Woods- 
Saxon profiles. The medium created in Pb+Pb colli- 
sions at Y^ijvw = 2.76 TeV is assumed to thcrmalize at 
to = 0.6 fm/c, at which time the temperature of the 
hottest point in central collisions is set to Tq = 520 McV, 
a factor of 1.3 larger than the highest temperature in 
Au-Au collisions at ^/sWn = 200 MeV at RHIC. This 
assumption is consistent with the measurement of the 
total multiplicity in Pb-|-Pb collisions at the LHC, which 
is a factor of 2.2 larger than that measured in Au+Au 
collisions at RHIC [2J]. The time evolution of medium 
is modeled by a one-dimensional boost-invariant cxpan- 
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FIG. 2: (Color online) The energy loss experienced by a 
100 GeV quark jet initially produced in the center of a head-on 
Pb+Pb collisions at ^/s^ = 2.76 TeV at the LHC. Contri- 
butions from various components are shown for comparison: 
total (solid line), vacuum radiation (dashed line), medium- 
induced radiation (dot-dashed line) and leading parton colli- 
sional energy loss (dotted line). 



sion, i.e., the temperature falls with time as t~^/^. The 
shower evolution is terminated when the medium temper- 
ature drops to Tc = 160 MeV, below which we assume 
the partonic shower suffers no energy loss in the medium. 

In Fig. [2] we show the medium- induced energy loss ex- 
perienced by a quark jet with initial energy E = 100 GeV 
produced in the center of the medium created in a head- 
on collision of two Pb nuclei at y^s]vjv = 2.76 GeV. We 
observe that, due to falling temperature and the varying 
transport coefficient q along the jet path, the length de- 
pendence of the energy loss become much weaker than for 
a static medium as shown in Fig. [TJ Such a jet, which is 
produced at the center of the collision region, experiences 
vacuum radiation before = to = 0.6 fm/c and then in- 
teracts with the medium for tf — Iq 5.5 fm/c before 
the medium temperature drops below Tc = 160 McV. 
The total energy loss from the jet cone during this pe- 
riod is AE « 22 GeV. Before t - to fa 1.7 fm/c, the 
energy loss is dominated by vacuum radiated gluons pro- 
duced before the medium evolution starts at to- After 
this time, medium-induced radiation takes over and fi- 
nally contributes about 70% of the total energy loss. 

In Fig. [21 we show our calculation of the medium mod- 
ification of the di-jet asymmetry factor Aj defined as 



A, 



E' 



T,l 



E' 



T,2 



Eta + Et,2 



(10) 



where ET,i,{i = 1,2) denotes the transverse energy of 
the leading and sub-leading jet, respectively. For back- 
to-back di-jet events in the vacuum, Aj is peaked at 
zero. The ATLAS Collaboration measured this quan- 
tity by requiring the trigger jet Et,i > 100 GeV and the 
second jet in the opposite hemisphere A(f) > 7r/2 with 
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FIG. 3: (Color online) Distribution of di-jet asymmetry factor 
Aj for p + p and Pb+Pb collisions at yijvjv ~ 2.76 TeV at 
the LHC. Left panel: 0-10% centrality; right panel: 10-20% 
centrality. 



Et,2 > 25 GeV. To proceed, we first generate vacuum 
di-jet events from PYTHIA [20| and obtain the distribu- 
tion for the di-jet asymmetry factor Aj in p + p events. 
The modification of each di-jet event in Pb+Pb collisions 
is obtained as follows. For each di-jet event, we sample 
its production points according to the distribution of the 
binary nucleon-nucleon collisions in collisions of two Pb 
nuclei. For asymmetric di-jets (Aj > 0.1), the trigger 
bias is taken into account by letting the higher energy jet 
propagate along the shorter path (implying a smaller en- 
ergy loss), and the other jet to propagate along the other 
direction. For nearly symmetric jet pairs {Aj < 0.1), 
such a trigger bias does not apply. 

As expected, the number of strongly asymmetric di- 
jets is significantly increased by the medium evolution 



which tends to let one jet lose more energy than the 
other due to the different path lengths of the two jets in 
the medium. The asymmetry of di-jets is more promi- 
nent in the most central Pb-fPb collisions (left panel 
of Fig. [3]) than in mid-central events (right). The de- 
pletion of energy inside the jet cone is a combination 
of collisional energy loss experienced by all shower par- 
tons, radiation outside the jet cone, and the scatter- 
ing of radiated gluons into angle outside the jet cone. 
From our fit to the data we obtain the average path- 
length weighted transport coefficient in central collisions 
{q) = {qL)/{L) = 1 GeV^/fm, where the average is over 
different production points and propagation directions. 
This corresponds to a value of g = 2.6 GeV^/fm at the 
highest temperature 400 MeV in Au+Au collisions at 
RHIC, consistent with the systematic analysis performed 
in Ref. m. 

In summary, we have studied the evolution of a jet 
shower propagating in a quark-gluon plasma and calcu- 
lated the loss of energy contained in a given cone an- 
gle. The medium modification of the shower spectrum 
and shape is described by a differential equation that 
incorporates both, collisional energy loss and transverse 
momentum broadening. Our approach provides a good 
description of the di-jet asymmetry observed by the AT- 
LAS Collaboration in Pb+Pb collisions at the LHC. The 
values of the parton transport coefficients are similar to 
those describing jet quenching at RHIC, extrapolated to 
the higher matter density at the LHC. This suggests that 
the quark-gluon plasma created at the LHC has similar 
properties as that studied by the RHIC experiments. 

This work was supported in part by Grants No. DE- 
FG02-05ER41367 and No. .de-scOOOSSQe, from the U.S. 
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